A RADIAL ANALOGUE OF POISSON'S SUMMATION FORMULA 



WITH APPLICATIONS TO POWDER DIFFRACTION 
AND PINWHEEL PATTERNS 
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Abstract. Diffraction images with continuous rotation symmetry arise from amorphous 
systems, but also from regular crystals when investigated by powder diffraction. On the 
theoretical side, pinwheel patterns and their higher dimensional generalisations display such 
symmetries as well, in spite of being perfectly ordered. We present first steps and results 
towards a general frame to investigate such systems, with emphasis on statistical properties 
that are helpful to understand and compare the diffraction images. An alternative sub- 
stitution rule for the pinwheel tiling, with two different prototiles, permits the derivation 
of several combinatorial and spectral properties of this still somewhat enigmatic example. 
These results are compared with properties of the square lattice and its powder diffraction. 



1. Introduction 

Since the discovery of quasicrystals some 20 years ago, mathematicians and physicists have 
gained a reasonable understanding of aperiodicahy ordered systems, in particular of those 
obtained from the projection method. Such sets are called cut and project sets, or model sets 
|17j . These are Delone sets of finite local complexity with respect to translations, which also 
means that any finite patch occurs in finitely many orientations only. 

Much less is known about aperiodically ordered systems with local patches occurring in 
infinitely many orientations, such as the pinwheel tiling of the plane or its three-dimensional 
counterpart, see [T^ and references therein. Arguably, these are closer to amorphous systems, 
but still perfectly ordered. To our knowledge, no mathematically satisfactory frame for the 
analysis of radially symmetric systems and the comparison of their spectral properties has 
been developed so far. It is the aim of this contribution to show first steps in this direction, 
by combining results and methods from discrete geometry with the more recent approach of 
mathematical diffraction theory. 

Our guiding examples are the square lattice and the pinwheel tiling of the plane. The 
diffraction of the square lattice is well understood, and it is not difficult to get some insight 
into its powder diffraction. The latter emerges from the presence of many grains in random 
mutual orientations, and thus requires a setting with circular symmetry. 

Circular symmetry is also a fundamental property of the pinwheel tiling, resp. its compact 
hull. This tiling has recently been reinvestigated from the autocorrelation and diffraction 
point of view ^Hl- However, hardly any explicit calculation exists in the literature, the reason 
being the enigmatic nature of the substitution generated pinwheel tiling. New insight is gained 
by means of an alternative construction on the basis of a substitution rule with two distinct 
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prototiles. This gives access to quantities such as frequencies (hence also to the frequency 
module), distance sets and the ring structure of the diffraction measure. 

In Section [21 we derive a radial analogue of Poisson's summation formula for tempered 
distributions, which is needed in the following analysis. In the last section, we outline the 
general structure by presenting a number of results. They are clearly distinguished according 
to their present status into theorems (with proofs or references), claims (with a sketch of 
the idea and a reference to future work) and observations (based on numerical or preliminary 
evidence). Since our alternative substitution will no doubt enable other developments as well, 
we hope that further progress is stimulated by the results presented in this paper. 

1.1. Notation and preliminaries. A rotation through q about the origin is denoted by 
-R^, and B^{x) is the closed ball of radius r centred in x. Whenever we speak of an absolute 
frequency (e.g., of a point set), we mean the average number per unit volume, whereas relative 
frequency of a subset of objects is used with respect to the entire number of objects. The 
absolute frequency of the points in a point set X (if it exists) is also called the density of 
X, denoted by dens(X). The set of non-negative integers is called Ng, and the unit circle 
is = {2 S C : \z\ = 1}. If M is a locally finite point set, the corresponding Dirac comb 
■~ SxeM'^x! where 6^ is the normalised point (or Dirac) measure in x, is a well-defined 
measure. The Fourier transform of g is denoted by g. Whenever we use Fourier transform for 
measures below, we are working in the framework of tempered distributions j2Uj . 

2. A RADIAL ANALOGUE OF POISSON'S SUMMATION FORMULA 

The diffraction pattern of an ideal crystal, supported on a point lattice F C M"^, can 
be obtained by the Poisson summation formula (PSF) for lattice measures or Dirac combs 

120113 El 

(1) 6r = dens{r)-6p*, 

where F* := {x G M*^ : x-y G Z for all y G F} is the dual lattice. The distribution- valued 
(or measure- valued) version follows from the classical PSF, compare [3 El? via applying 
it to a (compactly supported) Schwartz function. It is often used to derive the diffraction 
measure 7^ of a lattice periodic measure uj = g * 6p (with g a finite measure) . One obtains 
the autocorrelation measure 

(2) 7^ = dens(r) ■ig*g)* 6r, 

where g{g) := g{g) with complex conjugation 7 and g{x) = g{—x), and the diffraction measure 

% = (dens(r))2 • \g\^ ■ 6p*, 

see jl21 ^ for details. Our interest is to extend this approach to situations with circular (or 
spherical) symmetry. 

Let us first consider the square lattice Z^, and recall that its circular shells have radii 
precisely in the set 

(3) := {r > : r^ = m"^ + with m, n G Z} = {0, 1, ^/2, 2, \/5, 2^/2, 3, . . .}. 

This is the set of non-negative numbers whose squares are integers that contain primes p = 3 
mod 4 only to even powers. Moreover, on a shell of radius r G one finds finitely many 
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lattice points, their number being given by 



(4) ?7n(r) 



1, r = 0, 

4a(r2), r G \ {0}. 



Here, a(n) is the number of ideals of norm n in Z[i], the ring of Gaussian integers. This is 
a multiplicative arithmetic function, thus specified completely by its values at prime powers 
(see jni ^5 ) • They are given by 



1, ifp = 2, 

£ + 1, ifp = lmod4, 

0, if p = 3 mod 4 and i odd, 

1, if p = 3 mod 4 and i even. 



As is self-dual as a lattice, with dens(Z^) = 1, the PSF simplifies to 6^2 = (5^2 ■ 

Choose an irrational number a E (0, 1). By Weyl's lemma \TE\ , the sequence (na mod l)n>i 
is uniformly distributed in (0, 1). Consider the sequence (a;^)^>^ of measures defined by 

1 ^ 

(^) ^ ]V X] ^R"Z2 ' 

n=l 

where i? = R2na t'^^ rotation through 27rQ. If |x| = r > 0, the sequence is 
uniformly distributed on dB^(0), again by Weyl's lemma. Observe that all lattices R^l? 
share the same set of possible shell radii, namely of ©• This implies that, given an 
arbitrary compactly supported continuous function one has 

1 ^ 

^ 71=1 ^'^T^u 

where the measure is the normalised uniform distribution on dB^{0) = {x G : = r}, 
with fiQ = Sq. This establishes the following result. 

Proposition 1. The sequence (a;^)^>]^ of © converges in the vague topology, and 

a; := lim = V 7?n(r)^^, 

N^oo ^ — ' 

with shelling numbers ij^ (r) and probability measures /u^ as introduced above. □ 

It is obvious that the limit is, at the same time, also a limit of tempered distributions, i.e., 
a limit in 5'(M^). As the Fourier transform is continuous on 5'(M^), one has 

uj= { lim ujm] = lim Dm-. 

Kn^oo J N^oo 

Employing the ordinary PSF (pQ), one finds 

^ N ^ N ^ N 

n=l n=l n=l 
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where we used that (Rr)* = RF* for an isometry R, and dens(i?"Z2) = dens(Z2) = 1. 
Combining the last two equations, one sees that 

tD = Mm uJat = lim loat = uj. 

Using polar coordinates, the Fourier transform of /x^ — which is an analytic function 
because has compact support — can be expressed by a Bessel function of the first kind 
via the following calculation, 

Jl^ik)= e-2-^'=-^d;U,(x)= / — / e-2-i|'=l^?=°«^d^d(5,(^) 
iR2 Jo 27r Jq 

(6) 

Jo{27r\k\Q)d6r{g) = Jo{27r\k\r) 

where Jq{z) = X^^g (f)^^- This yields the identity 
(7) ^=Y1 ^□M^o(2vr|fc|r), 

to be understood in the distribution sense. It has the following consequence. 
Corollary 2. With r]^{r) and as introduced above, one has 

Yl ^□(^)^r = {Yl ^□('^)^r) = Yl ^oir)Jo{2TT\k\r), 
rexi^ rev^ rev^ 

where the last expression is to be understood in the distribution sense. □ 

Identities of this type can be viewed as measure-valued generalisations of classic Hardy- 
Landau- Voronoi summation formulae, see |131 Sec. 4.4] and references given there for details. 
However, in view of the rather delicate convergence properties, a direct verification via the 
PSF (0) seems a simpler approach. 

Observe that one can alternatively reduce the problem to one dimension and employ a 
Hankel transform, compare jl3l Sec. 4.4]. This requires a separate treatment of the trans- 
formed measure at the origin in fc-space. As this looks slightly artificial from the point of 
view of diffraction, we stick to ordinary Fourier transform. 

The version for a general lattice F CM.'^ reads as follows, where, in analogy to r]^{r), r}p{r) 
and r]p,[r) denote the number of points of F and F* on centred shells 9-6^(0) of radius r. 

Theorem 3 (Radial PSF). Let F be a lattice of full rank in , with dual lattice F* . If the 
sets of radii for non-empty shells are Dp and T>p,, with shelling numbers r]p{r) and r]p,{r), 
the classical PSF has the radial analogue 



(8) Vrir)l^r) = Y Vr{r)K =<iens{F) Y Vr'ir) l^r, 

where /i^ denotes the uniform probability measure on the sphere of radius r around the origin. 
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Figure 1. A plot of the radial structure of the function fii{k) of Q for 
dimensions d = \ (light grey), 2, 3, 4, and 5 (black). 



Proof. Select a sequence of isometrics (-R„)„>0) Rn ^ SO((i), such that {RnX)n>oi for a fixed 
X of length 1, is uniformly distributed on the unit sphere S'^. Consider then the sequence of 
measures defined by 



n=l 

The claim now follows from Weyl's lemma and the classical PSF in complete analogy to 
our previous planar example. □ 

The formula for general d can also be expressed in terms of Bessel functions of the first 
kind. Here, by standard calculations with spherical coordinates and integral representations 
of Bessel functions, one obtains 

r , d •^._i(2vr|A;|r) 

(9) Mk)= / e-2"*'=^d//^(x) = r(-' 



^ (7r|A:|r)i-i 



where T{x) is the gamma function and Ji_,{z) = {'^Y ^ £!rL+£+i') (f) • particular, 

1=0 ■ 



cos(27rA;r), if d = 1, 

Jo(27r|A:|r), if d = 2, 

sin(27r|fc|r) -f j _ q 

27r|fc|r ' II a - 



The analogue of Corollary |21 for d = 1 and F = 7L thus reads 

^ cos(27r/cm) = ^ = (5^, 

which is just another form of the ordinary PSF in this case (as radial averaging is trivial 
here). Figure ^ shows @ for r = 1 and various values of d. 

3. A SIMPLISTIC APPROACH TO POWDER DIFFRACTION 

The intensity distribution in powder diffraction emerges from a collection of grains in 
random and mutually uncorrelated orientations. Its precise theoretical description is difficult, 
compare |21| and references therein. 
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Here, we look into a rather simplistic approach that nevertheless captures the essence of 
the diffraction image. For simplicity, and for comparison with related pinwheel patterns, 
we explain this for the square lattice 1? . Instead of working with grains of finite size, we 
consider the superposition of entire lattices, with appropriate weights. Moreover, we make 
the restriction that there is a common rotation centre for all lattices, which we choose to be 
the origin. As a first step, let us take a look at 1? U Rl? ^ where R G S0(2) is a generic 
rotation (by which we mean that it is not an element of the group of coincidence rotations 
SOC(Z^) = S0(2,Q), see I3| for more on this concept). 

Lemma 4. Let R € S0(2) he a generic rotation, so that 1? n Rl? = {0}. Then, the 
autocorrelation of u = ^ {S^2 + Sj^i2 ] 



IS 



(10) Ico = + \^RI? + \\ 



where A is the Lebesgue measure in M . The diffraction measure of to is 

(11) %= E ^(^)'^- 

with 1(0) = 1 and I{x) = \ for all / x G U RI? . 

Proof. The first two terms in ()1U() are the autocorrelations of ^Z^ and of ^i?Z^. The third 
term originates from the cross-correlation between them, where we used the identity 

A = lim K-v^ 

r— >oo TTr 

yera2nB,.(o) 

with the limit being taken in the vague topology. This identity can be derived as follows. It is 
easy to see that each square in the square lattice is hit by x — y exactly once, if x is arbitrary 
but fixed and y runs through Rl? . (The exception is y = 0, but we can neglect it since it 
plays no role in the limit.) Thus each square in is hit with the same frequency in the 
limit. Since i? is a generic rotation, the sequence {x — y mod (1, 1)) is uniformly distributed 
in the fundamental domain [0, 1)^ of Z^, if x is arbitrary but fixed and y runs through Rl? . 
This is a consequence of the uniform distribution of {na mod l)„>i in (0, 1) if a is irrational. 
Thus, the points in the sum above are uniformly distributed in in the limit, and the series 
converges in the vague topology to the Lebesgue measure. 

Now, Equation (jllj) follows immediately from \ = and the PSF □ 

For completeness, let us briefiy comment on the situation that R E SOC(Z^). In this case, 
:= Z^ n RI? is a sublattice of Z^ of finite index, which is 1/ dens(6'). It is possible to show 
that w = ^ {bj2 + (5^22 ) has diffraction 

- -A 1a 1^ 

7a; - % + 4°Z2\0 + ^'^RZ2\e- 

The difference to the diffraction formula in Lemma |1] is concentrated on \ {0}, and hence 
plays no role in our further discussion when dens(6') — > 0, which happens under multiple 
coincidence intersections 0. We may thus disregard coincidence rotations for our present 
purposes. 
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Figure 2. Radial dependence of ring intensities of square lattice powder 
diffraction. The central intensity is not shown, see text for details. 
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Figure 3. Numerical approximation of Figure [21 based on summing Bessel 
functions, for radii r < 25. The required radial autocorrelation coefficients 
r]fj (r) are taken from Equation . 

Let us continue by considering multiple intersections. If all rotations are generic (in the 
sense that the lattices Ril? and Rjl? are distinct apart from the origin) and satisfy a uniform 
distribution property, we obtain the following result. 

Proposition 5. Let F = Tl? , F- = Rj'^'^ with generic Rj, i.e., F- H Fj^ = {0} for j ^ k, and 



define u = jj Ylj=i ■ Then, one has the identities 



7a. 



N - 1 
N 



^ 1 



N 



and 7^ 



iV- 1 1/1 v-^^ 



If {Rax)a^q is uniformly distributed on S"^ for some fixed x G we obtain 



1 ^ ^ 
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with %(r) as in @. 

Proof. The first statement — about — follows as in Lemma |^ Observing (Pj)* = RjT*, 
this implies 

where /(O) = 1 and I{x) = otherwise. Consequently, one finds 

TV- 1 1/1 \ 

By Theorem |31 in connection with Weyl's lemma and the self-duality of , the term in the 
brackets converges to X]r>o ('')^r as ^ oo. □ 

This simple argument shows that, after discarding the central intensity and multiplying the 
remainder by A^, one is left with a circular diffraction pattern in the spirit of the radial PSF 
in Theorem |31 Disregarding the central intensity, the shelling numbers ^/□(r) reflect the total 
intensity, integrated over the rings of radius r, in this idealized approach to powder diffraction. 
Consequently, in a measurement that displays the intensity along a given direction, the re- 
sulting radial dependence is given by -^2m^^ compare Figure 121 The numerical approximation 
in Figure 131 included for later comparison, shows a strong oscillatory behaviour, as a result 
of summing Bessel functions. Note that this approximation disregards positivity in favour of 
including circular symmetry. The comparison gives a good impression on the overshooting 
that originates from this approach. 

4. Application to tilings with statistical circular symmetry 

4.1. The pinwheel tiling. The prototile of the pinwheel tiling is a rectangular triangle with 
side length 1, 2 and \/^. The smallest angle in T is arctan(^). Here, we choose the prototile 
T with vertices (— ^, —^), (^) ~^)' f)' ™^ equip T with a control point at (0, 0). Every 
tile in a pinwheel tiling is either of the form R^T + x, or of the form R^ST + x for some 
X G M^, where S denotes the reflection in the horizontal axis and R^ rotation through a. The 
substitution a for the pinwheel tiling is shown in Figure (31 (left). By the action of cr, T is 
expanded, rotated by — arctan(i) and dissected into five triangles that are congruent to T. 
Formally, we define a{T) := {Qi{T) + xi, Q2{T) + X2, Q-i{T) + X3, Qi{T) + X4, Q^{T) + X5}, 
with Qi G 0(2) and xi G M^. The appropriate choices of Qi and Xi can be derived from the 
figure. In particular, <t(T) contains a triangle equal to T, thus we may choose Qi = 1 and 
xi = 0. The substitution a extends in a natural way to all isometric copies QT + x (where 
Q G 0(2) and x G M^) by a{QT + x) := Q(j{T) + Vbx. Since one of the triangles in a{T) 
equals T, the tiling PW can be defined as a fixed point of the substitution a. This results in 
a{PW) = PW = [Jn>o^"'i'^)- This convention follows |TH] but deviates from because 
it is advantageous for us to include the rotation in the definition of a. 

4.2. The kite domino tiling. It is advantageous to consider a second substitution which 
yields a closely related tiling. Consider the two prototiles K and D (for 'kite' and 'domino'). 
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Figure 4. The substitution for the pinwheel tihng (left) and for the kite 
domino tihng (right). The dashed hnes indicate how the dominos have to be 
dissected to obtain a pinwheel tiling. They are also needed to turn the global 
substitution rule for KD into a local one for its hull. 



where D is the rectangle with vertices {—\, ~\)^ {\} ~\)^ f)) f)) ^ quad- 
rangle with vertices (— ^), (^,— ^), (^,^), (~"^)|)- Both tiles consist of two copies of 
the pinwheel triangle T, glued together along their long edges. Every pinwheel tiling gives 
rise to a kite domino tiling by deleting all edges of length The substitution for the kite 
domino tiling is shown in Figure |3 (right). Essentially, it is o"^ applied to the new prototiles, 
where a is the substitution for the pinwheel tiling. If q denotes this new substitution, a fixed 
point is given by KD = IJn>o B^iD)- This is a global substitution that works on the specific 
D defined above, compare [H] for a discussion of the general substitution concept. The two 
tilings, PW and KD, are mutually locally derivable (MLD) in the sense of 0. Essentially, 
this means that one is obtained from the other by local replacement rules and vice versa. 
These rules are evident from Figure [l] 

As PW and KD are MLD, the global substitution q also induces a local substitution in 
the sense of p], where two (oriented) copies of the domino have to be distinguished. Loosely 
speaking, this follows from the observation that the local surrounding of any domino in KD 
or one of the other elements in the hull defined by KD determines the type of the domino, 
and hence how to apply the substitution to it. 

We equip K with control points at (0, 0), (|, -g), and D with control points at (0, 0), (0, 1). 
Then, the set of all control points in KD is equal to the set of control points in PW . This 
specific discrete point set is a Delone set, denoted by A in the sequel. It is not hard to see 
that A is MLD with both PW and KD, even though one direction of the replacement rule is 
less obvious than the one linking PW and KD. 

Observe that the relative and absolute frequencies of triangles T in the pinwheel tiling are 
both equal to 1. For the relative frequencies, this is clear since there is only one prototile. The 
absolute frequency then follows from the fact that this prototile has area 1. Every triangle 
of PW carries exactly one control point, so dens(yl) = 1, too. For the kite domino tiling, 
standard Perron- Frobenius theory yields the relative frequencies of kite (resp. domino) as 
^ (resp. j^). Thus, the absolute frequency of a kite (resp. domino) is ^ (resp. ^). The 
substitution matrix can be extracted from Figure |^ 

Discrete structures that are MLD lead to dynamical systems that are topologically con- 
jugate |14) . Therefore, we formulate the hull on the basis of the Delone sets of the control 
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Figure 5. A patch from a kite domino tiling. The control points of the white 
tiles are contained in H? . The control points in the other tiles are contained 
in rotated copies of Z^, where the common rotation centre is indicated in the 
figure. The possible rotations are described in Claim ^ 



points. Let A be one such set, e.g., the set A of control points of PW as defined above. Define 
the orbit closure 

I p-r 

X(yl) =R2+/^' 

in the local rubber topology (LRT) 4,, which is compact. This topology is slightly different 
from the one introduced in |19|. But in the present case, both topologies yield the same hull 
and are equivalent on X(yl). 

The following result is well-known jl9l I18j . 

Theorem 6. The hull X(yl) is -symmetric. Moreover, {'K[A),'M?) is a strictly ergodic 
dynamical system, i.e., it is uniquely ergodic and minimal. All elements of X(yl) possess the 
same autocorrelation measure J = Jy[ and the same diffraction measure 7^. Both measures 
are -symmetric. □ 

In this sense, speaking of the diffraction of the pinwheel tiling or its control points has a 
unique meaning. At this point, we know that 

7 = (dens(yl))^ (5o + (7)cont = -^o + (7)cont 

because a translation bounded circularly symmetric measure cannot contain Bragg peaks 
other than the trivial one at /c = 0. That the intensity coefficient of 5q is (dens(yl))^, hence 1 
in our case, is a consequence of the equation 

7({0}) = lim -^^{BM) = ( lim ^ <5^(i?.(0)) V = {dens{A)f . 

It can be proved by means of a Fourier concentration argument in connection with results 
from alternatively, see ^21 for another approach and JH] for details on the case at hand. 
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It remains to determine the structure of (7)i-ont ™ore closely, in particular the separation into 
singular and absolutely continuous components. 

4.3. Results and observations. The following claims and conjectures are stated for the 
particular tiling KD or for its set of control points A. Since mutual local derivability extends 
to entire hulls, we remain in the situation of Theorem |HJ Proofs of the claims are either 
outlined here or will appear in 2\ Let us start with our main observation. 

Observation 1. The diffraction measure 7 of the pinwheel tiling is of the form 

1 = ^0 + (7)sc + (7)ac- 

There is a countable set T>* C M>o such that 

(7)sc = I{r)l^r- 
r&V*\{0} 

The set V* seems to he locally finite, i.e., discrete and closed. Moreover, (7)ac seems to he 
non-vanishing. 

The set V* is a subset of another set, T> = P^, which shows up in the determination of the 
autocorrelation and is described in more detail below. 

Claim 1. Let A be the set of all control points of KD. Then, A is a Delone subset of a 
countable union of rotated square lattices. More precisely, with 6 = 2arctan(^), 

Aci\J R^I?). 

Using the geometry of the kite domino tiling, it is not difficult to establish this property. 
Figure |S1 may serve as a visualisation on a small scale. This structure suggests that the 
diffraction of the pinwheel tiling might share some features with that of the powder diffraction 
of 1?, see the comment after the proof of Lemma Since Rq acts as multiplication by 
id ) ' next result is immediate. Note that Rg is well known from the coincidence site 
lattice problem of the square lattice, compare [Sj, and has played a crucial rule in the explicit 
calculations in jl8j . 

Claim 2. Let A he the set of all control points of KD. It satisfies 

{(f ,f ) :m,nEZ, G Nq}, 



and the distance set V = T>y^ := {\x — y\ : x,y £ A} is a subset of {y : p,q,£ & Nq}. It 

is the same set for all elements of X(yl). □ 

We believe that this subset relation is quite sharp, i.e., that the difference between the two 
sets is not too large, in the following sense. 

Claim 3. All values r = \/ p"^ + q"^ with p, Q G Ng occur in T>. Moreover, for each £ £ Nq, 
there are infinitely many p,q such that \J ^ ^'^ is contained in T>. 

In fact, numerical computations suggest the following, stronger property. 



Observation 2. For each £ G Nq, T) contains all but finitely many values of the form y ^^-^ 
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In the sequel, the absolute frequencies of distances are of interest, wherefore we define the 
radial autocorrelation function 

1 



(12) 



r]{r) 



lim 

i?->oo 



\x-y\=r 

The limit exists due to unique ergodicity, see Theorem El 
autocorrelation of yl — and hence of X(yl) — as 



1. 



This permits us to write the 



(13) 



7 



re© 



r)lJ'r 



with fj,j. as in Section [2j This follows from Theorem El together with Claim El Note that A is 
repetitive ^Hl (defined up to congruence), so that all r](r) in (fTTH) are strictly positive, by the 
definition of P in Claim [2 Clearly, 7 is both a positive and positive definite measure, and a 
tempered distribution on M?. 

The following theorem holds in more generality than the context of this paper. In fact, it 
holds for all substitution tilings which are of finite local complexity (FLC) and self-similar. 
The latter means that AT- = UTeo-(T.)^- Roughly speaking, this means that the support 
of the substitution of each prototile T- is similar to T-. Finite local complexity means that, 
for some i? > 0, the tiling contains only finitely many local patches of diameter less than R, 
up to congruence. Note that FLC is frequently defined with respect to translations, whereas 
we define FLC here with respect to congruence, because of the nature of the pinwheel tiling. 
Both properties, FLC and self-similarity, hold for the majority of substitution tilings in the 
literature IIUJ. 

Theorem 7. In any self-similar substitution tiling of finite local complexity [w.r.t. congru- 
ence) with substitution factor \, all relative frequencies are contained in Q(A'^). Moreover, 



the tiling can be scaled such that all absolute frequencies are contained in Q(A'') as well. In 
particular, if X'^ is an integer, and if the tiling is appropriately scaled, all relative and absolute 
frequencies are rational. 

For a proof, we refer to [2] . This theorem, applied to the pinwheel tiling or the kite domino 
tiling, yields that all frequencies are rational. 

Claim 4. All values 0/ < 5 together with rj{r) are given in the following table. Values 
marked with an asterisk are numerically based conjectures, all other values are exact. 



r2 





1 

5 


1 


8 
5 


9 
5 


49 
25 


2 


13 

5 


81 
25 


17 
5 


4 


113 
25 


5 


r]{r) 


1 


5 
11 


439 

165 


1 
2 


67 
165 


4 
165 


7* 
2 


142 
165 


4 
165 


10* 
11 


3* 


8 * 
165 


73* 
15 



The value 77(0) is the absolute frequency of the control points. Since each triangle in the 
pinwheel tiling has area 1, we have r/(0) = 1. The distance |x — y| = l/\/5 occurs precisely 
once within each kite. Kites have absolute frequency but the distance — y| has to be 
counted twice in view of formula ((T^ . hence r]{l/^/5) = ^. The other exact values contained 
in the table can be established similarly, but require more sophistication. With some further 
effort, one can determine the frequency module of A, which is the Z-span of the absolute 
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1 6 1 6 1 2 5 17 3 2 1 



5 55 22 55 22 33 264 220 110 165 264 

Figure 6. The 11 vertex stars of the pinwheel tiling (up to congruence) and 
their absolute frequencies (i.e., frequencies per unit area). 

frequencies of all finite subsets of A, the latter standardised to having density 1 (which is our 
natural setting here). 

Claim 5. The frequency module of A, and hence of X(yl), is J- = { 264" 5^ | ?7i G Z, ^ G Ng}. 
In particular, it is countahly, though not finitely generated. 

This result is derived from the absolute frequencies of the vertex stars in the pinwheel 
tiling, see Figure El These, in turn, can be derived from the frequencies of kites and darts in 
the kite domino tiling. Details will be given in 

We proceed by considering the Fourier transform in relation to the distance set V. Since the 
Fourier transform is continuous on the space 5'(M^) of tempered distributions, H13|l becomes 

(14) = Vir) IJ-r) = Yl ^^^^ 

where the sum is to be understood in the distribution sense. By Bochner's theorem, 7^ is 
again a positive and positive definite measure, and the equation can be understood as a vague 
limit as well. 

Let us apply Theorem |31 to the pinwheel pattern A. According to Claim the distances 
in P arise from lattices of the form ^^Z^ with £ E Nq. The corresponding dual lattices, 
which are relevant for the diffraction analysis, are h^l"^!?^ which have distance sets 
Since Tl? contains a lattice congruent to ^/5Z2, it follows that h^l'^V^ C for all £ G Nq, 
and that T)^ contains all distances that seem relevant from this point of view. At this stage, 
we have not found a compelling argument to exclude the possibility that certain subsets of A 
contribute relevant distances in a coherent fashion, and thus - by unique ergodicity - further 
rings I(r)//r (with r P^) to 7. However, on the basis of Theorem|3J it is at least plausible 
that the set P* is closely related with P^, or even equal to it. In this case, 7^ would show 
rings /(r) /i^, for r G P^ only. According to Claim |31 all distances have positive frequencies, 
so that this assumption would lead to 

(7)sing = Y ^ir) fJ'r = ^0 + (7)sc = ^0 + ^^^^ 

with 1(0) = 1 and I{r) > for ah r E P^ . 

Let us continue our discussion on the basis of this hypothesis. If ObservationElholds, then, 

for each £ > 0, only finitely many distances of type v/^~^^ ^^'^ missing, each one contributing 
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m 








Figure 7. Numerical approximation of the radial intensity dependence of 
the pinwheel diffraction pattern. It is based on Equations © and (|14|) . dis- 
regarding contributions to the central intensity. The radial autocorrelation 
coefficients r}{r) are estimated from a patch of radius of about 56 (thus ef- 
fectively using the fifth iteration of the substitution a). The vertical scale 
is arbitrary in the sense that it has no meaning without local averaging and 
integration. 



Figure 8. Detailed view of the shoulder to the left of the peak at A; = 1, 
calculated from the sixth iteration of the substitution rule, with the large peak 
intensity truncated. 

an absolutely continuous part to 7. Since 7 exists as a translation bounded measure, and since 
75i^g is already covered by the results above, every additional contribution has to contribute 
to 7ac- This motivates the conjecture that (7)30 7^ 0. 

Observation n is also supported by numerical computations of 7. A naive numerical anal- 
ysis of a large finite portion of A would yield no relevant result due to the nature of the 
pinwheel pattern. In particular, the number of orientations grows only logarithmically with 
the radius, while the number of tiles grows quadratically. The results above allow a more 
meaningful computation. In particular, by employing the S^-symmetry, the problem becomes 
one-dimensional, and by ((Tj) , the Fourier transform is expressed explicitly as a sum of weighted 
Bessel functions. As before, this approach disregards positivity of the intensity function, and 
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the result displayed in Figure [7| shows strong oscillations with significant overshooting, similar 
in kind to the ones observed in Figure El No smoothing of any kind was used. 

A comparison of the two diffraction images (Figures |31 versus E|) also supports our claim 
about the possible radii of pinwheel diffraction rings. One noticeable difference is the higher 
peak at A; = \/5, which is due to the pairs of points in A with distance 1/^/5. The existence 
of positive shoulders in Figure [7| (such as that to the left of the first peak at /c = 1 , a blow-up 
of which is shown in Figure EJ is another significant difference to Figure El and is one of 
the reasons why we expect a non- vanishing radially continuous contribution, hence giving an 
absolutely continuous component to 7. 
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